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Abstract 



The modified Gross-Pitaevskii equation was derived and solved to obtain the ID solution 
in the zero-energy limit. This stationary solution could account for the dominated contri- 
butions due to the kinetic effect as well as the chemical potential in inhomogeneous Bose 
gases. 
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Studies of collision phenomena in rather cold gases, e.g., dilute Bose gases, have recently at- 
tracted many researchers' attention [1-2]. One relevant research interest is about the solution 
of the appropriate and modified Gross-Pitaevskii equations for different dimensions [3-4]. New 
possibilities for observation of macroscopic quantum phenomena arises because of the recent re- 
alization of Bose-Einstein condensation in atomic gases [1-2]. There are two important features 
of the system - weak interaction and significant spatial inhomogeneity. Because of this inhomo- 
geneity a non-trivial zeroth-order theory exists, compared to the first-order Bogoliubov theory. 
This theory is based on the mean-field Gross-Pitaevskii equation for the condensate if) -function. 
The equation is classical in its essence but contains the (h/2p) constant explicitly. Phenomena 
such as collective modes, interference, tunneling, Josephson-like current and quantized vortex 
lines can be described using this equation. The study of deviations from the zeroth-order theory 
arising from zero-point and thermal fluctuations is also of great interest [5-7]. Thermal fluctua- 
tions are described by elementary excitations which define the thermodynamic behaviour of the 
system and result in Landau-type damping of collective modes. 

As a preliminary attempt, followling the mean- field approximate formulation in [3], in this let- 
ter, we plan to investigate the ID solution for the modified Gross-Pitaevski equation in the 
zero-energy limit. This presentation will give more clues to the studies of the quantum non- 
equilibrium thermodynamics in inhomogeneous (dilute) Bose gases and the possible appearance 
of the kinetic mechanism before and/or after Bose-Einstein condensation which is directly linked 
to the particles (number) density and their energy states or chemical potentials. 
The generalization of the Bogoliubov prescription [8] for the ^-operator to the case of a spatially 
nonuniform system is 





where ipo is the condensate wave function. This is an expression of the second quantization 
(^-operator) for atoms as no = Nq/V = \ipo\ 2 and (j) ^ \/no, Nq is the number of atoms in 
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the condensate. To neglect <p means neglecting all correlations and this is a poor approximation 
when distances between particles are of the order of the effective radius r e of the atom-atom 
interaction. To overcome this problem the assumption that the atomic gas is dilute : n r\ <C 1 
should be made [1-3] (n = N/V, N is the number of atoms confined in V). We can use the 
procedure of the quantum virial expansion to calculate the energy of the system. We have the 
form 

E = E, + 9 -j n 2 (r)dr, (2) 

for the energy of slow particles, where Eq is the energy of the gas without interaction, n(r) os 
the density of the gas, g = 4-irah 2 /m, a is the s-wave scattering length and m is the particle's 
mass [3]. After taking into account the correlations (in above equation so that we can neglect 
4>) and considering E as an effective Hamiltonian we then have the celebrated Gross-Pitaevskii 
equation 

d h 2 V 2 

ih-Mr,t) = [--^f- + U) + g\Mr,t)\ 2 ]M',t) (3) 

which describe the dynamics of a non-uniform non-ideal Bose gas at T = 0. Here, V ex t(r) is 
the confining potential. If the gas is in its ground state, the time dependence of V'o is given by 
ipo ~ exp(—ifj,t/h), where \i is the chemical potential of the gas [1-3]. We thus obtain the form 

[ ~^Sr + Vext{r) + ^ o(r ' t)|2 " ^ o(r ' t] = (4) 

which could be stationary once t is fixed or selected. We noticed that an equation of above form 
has been considered before in connection with the theory of superfluidity of liquid helium close 
to the A-point [9]. 

To investigate our interest here, we shall consider the ID solution of equation (4) for the case 
of zero-energy limit. Firstly we consider a 3D Bose gas confined tightly in one dimension and 
weakly in the remaining two dimensions on a length scale It {=\Jfi/2muj for a harmonic trap of 
angular frequency of a;). A collision between two condensate particles will typically occur over 
the characteristic length scale l co i once k is much larger than l co i (thus we can use a local density 
approximation) . 

We now model the pair wavefunction of two atoms in the medium by that a single particle with 
the reduced mass moving in a potential which consists of a circularly symmetric box of radius 
R and a hard sphere of radius R a located in the centre of the box. Following the reasoning in 
the derivation of equation (1), we introduce a similar bias or ghost-effect for (f> : * which can 
be relevant to certain critical or kinetic (non-equilibrium) effect (or configurational dissipation) 
[1-2,5,7,10-11] so that fi^f = e in the zero-energy, zero-momentum limit. It is presumed that 
e still reaches zero in the homogeneous limit. The problem for a 2D tp(r,9) becomes, after 
referencing to the bias 

d 2 d d 2 
dr 2 rdr d9 2 
and, in fact, for a circularly symmetric tp( r )i 

r d 2 d , , 
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with the boundary conditions : the wave function vanishes on the inner radius (tp = as r = R a ), 
and reaches an asymptotic value at the edge of the box (i/j — > IT as r = R). 
We can obtain the solution 

where i? a < r < R. The extra energy caused by the curvature of this wave function resulting 
from the scattering potential is 

h 2 ,2tt rR 2 h 2 rr R 3 -Rf, R 2 -R 2 n/ l R-R a „ 

AE= / / \Vtb(r ) 2 rdrd6 = Tr— \e[ ^ + e - ° ( r-±) + 

2mJo J Ra l ^ y n m x 1 12 ln(R/R a ) V 16 4 ;J 

;i [ n + £ (i?- J R a )(^-^-i)]}. (8) 
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Note that, this energy depends upon the size of the box R, which is indeed the length scale 
relevant for the scattering of two particles in two dimensions [1-3]. The scattering of two particles 
in a many-body system should obviously not depend on the size of the system as a whole when 
R becomes large, and so we must interpret R as the physical relevant length scale l co \. The 
appropriate length scale over which a many-body wavefunction changes is the healing length 1^, 



given in homogeneous Bose condensed systems by lh = h/ \f2m §2D no = ft,/ \f2mjl [1-3], and so 
it is this which must be used in equation (8). We recall that |I1| 2 corresponds to the condensate 
density no and the homogeneous limit for a pair interaction strength can still be recovered from 
above equation. 

If we plot tp(r) w.r.t r (in terms of units of l t ) for different e (0.0005, 0.001, 0.005, 0.01) and 
the same II (0.01) into a figure then the subsequent presentation shows the significant effect of 
the kinetic part (say, II) due to e. The effect of boundary conditions, like II is minor. From the 
definition of e = u^, we can understand that the contribution of the chemical potential for ip 
of the inhomogeneous gases is indeed dominated, too. We shall investigate more complicated 
problems [12-13] by the same approach in the future. 
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